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We report on surface wave pattern formation in a Fara- 
day experiment operated at a very shallow filling level, where 
modes with a subharmonic and harmonic time dependence 
interact. Associated with this distinct temporal behavior are 
different pattern selection mechanisms, favoring squares or 
hexagons, respectively. In a series of bifurcations running 
through a pair of superlattices the surface wave pattern trans- 
forms between the two incompatible symmetries. The close 
analogy to 2D and 3D crystallography is pointed out. 

PACS: 47.54.+r 47.20.Ma 47.20.Lz 



When a fluid layer is vibrated vertically patterns of 
standing waves occur at the liquid-air interface. This 
experiment, first studied by Faraday Q in 1831, has be- 
come a paradigm for the investigation of spontaneous 
pattern formation From an experimental point of 
view, the Faraday setup is particularly attractive since 
the characteristic length and time scales are under ex- 
ternal control and relaxation times are pretty short. If 
the vibration signal is sinusoidal and weakly supercriti- 
cal one generically observes well ordered standing wave 
patterns in the form of squares or lines, oscillating with 
half the frequency of the excitation (subharmonic re- 
sponse). By carefully tuning the drive frequency within 
the crossover regime between gravity and capillary sur- 
face waves, more complicated (quasi-periodic) patterns 
with an 8 or even 10-fold point symmetry have been pre- 
dicted and observed Edwards and Fauve intro- 
duced the idea of a two-frequency excitation signal. The 
resulting interplay between the competing modes gave 
rise to unexpected new phenomena like hexagonal, trian- 
gular, quasi-periodic, spatially localized structures and 
even superlattice type patterns ||,^. However, the large 
number of control parameters of a multi-frequency drive 
signal renders such an experiment problematic as far as 
a systematic exploration of parameter space is concerned 
and it lacks an intuitive understanding of the fundamen- 
tal principles. By use of a viscoelastic liquid, it has been 
shown 1^ that the familiar subharmonic Faraday reso- 
nance can be preempted by a synchronous (harmonic) 
response, leading to a bicriticality as well. However, 
such fluids are more complicated than Newtonian liq- 
uids and the characterization of their nonlinear rheology 
is incomplete. For Newtonian fluids the parameter re- 
gion in which the response switches from subharmonic 
to harmonic is difficult to access |^-|ll[| since it requires 



a very shallow filling depth, large shaking elevations and 
thus a powerful vibrator. This article describes a system- 
atic investigation of a Newtonian liquid under a single- 
frequency drive operated close to the bicritical crossover. 
The interesting new feature is a "phase transition" from 
a quadratic pattern to a structure of hexagonal symme- 
try. The symmetry change proceeds via a sequence of 
superlattices with 2-fold and 6-fold point symmetry. 

Our new shaker system has a force of 4800 N and a 
maximum peak-peak elevation of 5.4 cm. It is operated 
at a rather low drive frequency of 8— lOi^z. The combina- 
tion of these parameters with a filling depth of 0.7mm, 
only, is necessary to access the harmonic-subharmonic 
bicriticality. The container with an inner diameter of 
290mm (equivalent to about 15 times the wavelengths) 
is sealed by a glass plate and temperature controlled at 
25° ± 0.1°C. At this temperature the sample fluid (low 
viscosity Silicon oil, Dow Corning 200) is specified by a 
kinematic viscosity v = 9.7 x 10~^m^/s, a surface ten- 
sion of <j = 0.0201A/'/m and a density p = 934/c^/m^. 
The actual acceleration experienced by the container is 
recorded by a piezoelectric device. A feedback loop con- 
trol limits disturbing anharmonicities to a level of less 
than 0.2%. Our visualization method for the evaluation 
of pattern symmetries has been described elsewhere [Q. 
Although this simple light reflection technique provides 
high contrast pictures, the relation between the recorded 
intensity and the surface profile is not trivial. Neverthe- 
less for simple structures like squares, we were able to 
reconstruct the profile by the following procedure: Start- 
ing from an estimated surface profile composed of a small 
number of spatial Fourier modes, we computed the light 
distribution of the expected video image by means of a 
ray tracing algorithm. Then we adapted the mode ampli- 
tudes and their relative phases to optimize the agreement 
between the calculated and recorded video pictures. 

The linear stability theory Q evaluated for our fluid 
predicts that the bicritical threshold, at which the sys- 
tem changes from the harmonic instability at lower to 
the subharmonic instability at higher frequencies, occurs 
at a drive frequency of VLb/'^'t^ = S.THz. This is in 
good agreement with the present results. Fig. |l| depicts 
a calculated neutral stability diagram (acceleration am- 
plitude a vs. wavenumber k) for an excitation frequency 
Q/27T = 9.5Hz slightly above Qb- The minimum of the 
left (right) resonance tongue defines the wavenumber kg 
(kh) of the subharmonic S (harmonic H) instability. The 
absolute minimum ac corresponds to the onset of Fara- 
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day waves. Fig. || presents a phase diagram classifying 
the surface patterns as observed during a quasistatic am- 
phtude ramp at fixed drive frequencies. Starting from 
a subcritical drive e = a/ac — 1 = —2% the amplitude 
a is increased by steps of 0.2%. After each increment 
the scan is suspended for 240 sec and a surface picture is 
taken thereafter. Having reached e = 10% the ramp is re- 
versed. There is no noticeable hysteresis for the primary 
onset between upward and downward scans. 

Within the drive frequency region of the harmonic 
Faraday instability {Q < Qb) the bifurcation scenario 
is rather conventional: Entering subregion V from below 
(see Fig. ||), we find a perfect hexagonal surface tiling 
which persists up to the maximum drive amplitude. 

The focus of the present paper is on the subharmonic 
region Q > Qb^ where the onset pattern (region //) is 
quadratic, while hexagons (region V) occur at rather el- 
evated shaking amplitudes as a higher bifurcation. Our 
aim is to describe the bifurcation sequence, which results 
from an amplitude scan at the fixed drive frequency of 
Q = 9.5 Hz (gravity wave regime). The primary Fara- 
day pattern (region // in Fig. ||) has a subharmonic 
time dependence and exhibits a perfect quadratic sym- 
metry as shown in Fig. ||a. The associated spatial power 
spectrum (Fig. ||b) indicates the fundamental wave vec- 
tors k^i and k5'2 but also pronounced contributions from 
higher harmonics, in particular k^i + k5'2- The appear- 
ance of square patterns in low viscosity gravity waves 
agrees with a small amplitude theory expanded around 
the subharmonic instability threshold . At those small 
values of e the competing harmonic Faraday modes do 
not noticeably affect the pattern selection process. The 
shallow filling level used in our setup makes the surface 
elevation profile very anharmonic with high, narrow tips 
but broad, shallow hollows. We have re-constructed the 
surface profile belonging to Fig. ||a by our ray tracing 
algorithm. The spatial dependence of the interface de- 
formation is decomposed according to 

7^(r) = ^AiCos{ki -r + ^i), (1) 

i 

where k^ G {k^i, ks2, (k^i ± ^32), 2ksi, 2ks2, 2{ksi ± 
ks'2)} and r = {x^y). A density plot with the gray 
level proportional to the local surface elevation is given 
in Fig. |3|d together with the related video image (c) as 
computed by ray tracing the refiected light. 

With increasing drive amplitude the harmonic Fara- 
day instability gradually gains infiuence upon the pat- 
tern selection dynamics. When entering region /// we 
observe a continuous transition to a \/2 x superlat- 
tice (see Fig. ^,b), characterized by a new subharmonic 
mode with wave vector \<.d1' Even though the associated 
Fourier spectrum also shows a second mode ]<.d2 orthog- 
onal to the first one, the amplitude ratio Adi/Ab2 — 4 
indicates that k/)i is strongly prevailing. Therefore the 
original quadratic invariance of the pattern is broken and 



replaced by the simpler rectangular symmetry. The sur- 
face profile as derived by ray tracing suggests that the 
mode liDi enters Eq. |l| in the form cos (k/^)! • r — 7r/2). 
Assuming (j)si = (/>s'2 = (by a proper choice of the ori- 
gin) the cubic nonlinearity AsiAs2A'^-^ enters the am- 
plitude equation for Adi and thus restricts the relative 
phase (j)Di to or 7r/2, depending on the sign of the 
related coupling coefficient. The former (latter) value 
is associated with a modulative (displacive) mode. De- 
tails will be published elsewhere. For our system the 
second case applies, which becomes evident by compar- 
ing the video images Figs. |a and |a: With mcreasms 
order parameter Adi the rows of elevation maxima (con- 
nected by dashed lines in Fig. ||a) are displaced in oppo- 
site directions as indicated by the arrows. The displacive 
mode ]<.Di is excited by a pair of wavevector triads ac- 
cording to the geometrical relations = — ^si 
and ]<.Di = ^H2 — ^S2' The modes I^hi^ ^H2 oscil- 
late synchronously with the drive. Their contribution to 
the triad-couplings becomes energetically efficient since 
\^Hi\ = |ki^2l — \/5/2 :^ 1.58 coincides almost exactly 
with the wavenumber kh — 1.59 associated with the har- 
monic Faraday instability (Fig.l). 

The transition from region /// to IV is accompa- 
nied by a very slow rearrangement of the pattern and 
the occurrence of defects. A snapshot taken during 
this crossover process is shown in Fig. As e is in- 
creased, the displacive mode ]<.di gradually dies out and 
thus liHi and ]<.h2 enter into a new triad together with 
the nonlinear harmonic mode k/^s = k5'2 — k^i. Even 
though l^Hsl = a/2 does not exactly match the unstable 
wavenumber band around kh (see Fig. ^ the drive am- 
plitude is apparently high enough to allow this detuning. 
The six-armed star of wavevectors 'kni shown in Fig. ^ 
already suggests a 6-fold rotational invariance, but the 
true symmetry of this " quasi- hexagonal" transient state 
is still 2-fold: The angles between ]<.hi and ]<.h2 and be- 
tween ]<.H2 and kji/a are 56° and 62°, respectively, rather 
than 60°. /^From the horizontal streaks of the Fourier 
spectrum of Fig. ^ it can be seen that the translational 
coherence in x-direction is in a process of disintegration, 
while the spatial periodicity in ^/-direction and the 2-fold 
point symmetry are still preserved. Nevertheless the mu- 
tual resonance among the ]<.Hi is the precursor of the 
hexagonal symmetry, which will follow. 

The transient re-orientation process comes to a halt 
when the pattern has accomplished the ideal hexago- 
nal symmetry (region IV and Figs. ^,f). At this stage 
the harmonic Faraday modes govern the pattern selec- 
tion process. The structure depicted in Figs. ^,f can be 
understood as a >/3 x >/3 superlattice of the hexagonal 
lattice of the final state (region V) . Apparently, the sub- 
harmonic contributions to the Fourier spectrum are dy- 
namically slaved, since the 6-fold symmetry dictated by 
the ki:/-modes is recovered in number as well as in ori- 
entation by the new subharmonic set {k^i, k5'2, ^ss}- 
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Note that the subharmonic time dependence does not 
ahow a mutual resonance among the k^^. The relative 
30°-orientation between the S-star and the H-star results 
from the triad cross-coupling between harmonic and sub- 
harmonic modes. The appropriate geometrical resonance 
conditions k^i +k5'2 = ^hi^ etc., also enforce the length 
of the subharmonic wavevectors to reduce from unity to 
\^Si\ = \/5/6 :^ 0.91. It can be shown that the funda- 
mental modes contributing to the x superlattice 
enter Eq. ^with equal spatial phases (j)Hi = (j^Si = 0. The 
equality (j)Hi = (I>H2 = (t>H3 = follows from the mutual 
resonance among the kji^-modes, while the cross-coupling 
between k^-- and kji^-modes enforces <psi = 052 = 4^33- 
The remaining condition (psi = is more subtle and es- 
tablished only at quintic order in the associated ampli- 
tude equation. With the reasoning outlined in Ref. Q 
one obtains either + (j)s2 + 053 = or 7r/2. Since the 
former (latter) case is associated with a 6-fold (3-fold) 
symmetry, we conclude from Fig.||e that (j)si = applies 
in our system. 

If the drive amplitude is rised into phase region V the 
long wavelength modulation of the x ^/3 superlattice 
disappears and the pure hexagonal state as depicted in 
Fig. 1^ survives. Simultaneously the r^/2-component in 
the temporal power spectrum of the surface oscillation 
dies out. 

To summarize: The present letter reports a Faraday 
experiment in a very thin fluid layer. Drive frequency 
and filling depth are adjusted such that waves with sub- 
harmonic and harmonic time dependence become simul- 
taneously unstable. The different time dependencies im- 
ply distinct wavelengths, but they are also responsible for 
different nonlinear pattern selection mechanisms, which 
favor either squares or hexagons. In a series of bifurca- 
tions the transition from one symmetry to the other takes 
place in a surprisingly coherent manner running through 
a displacive \/2 x and a v^S x ^/3 superlattice. 

A comparison to crystallographic phase transitions is 
in order. Hexagonal structures on one side and square 
(in 2D) or cubic (in 3D) structures on the other side 
are incompatible since there is no group-subgroup rela- 
tion connecting the different space groups. Hence such 
transitions involve a reconstruction of the lattice via the 
formation of lattice defects. The hcp-fcc transition, oc- 
curing e.g. in solid 4He, is a prominent 3D example. 
Here the relevant defects are stacking faults. The same 
principles hold of course for the present patterns. At the 
same stage along the phase sequence from // to IV there 
has to be a reconstructive transition involving defects. In 
fact the pattern of Fig. ||c taken just at the /// — IV- 
boundary shows such as line defect (running in the verti- 
cal direction in the right half of the picture). In contrast 
to the hcp-fcc transition, the transition between the two 
incompatible symmetries of our system does not occur in 
a single step, but involves two intermediate phases. Tran- 



sition of the type // — /// and V — IV are also known in 
2D crystallography. The transition from a simple hexag- 
onal lattice to a ^/3 x ^/3 superstructure has been ob- 
served for instance in monolayers of adsorbed on 
graphite [l4j . The displacive transition // — /// is analo- 
gous to the reconstruction of the (100) surface of Wolfram 
crystals Ref. [|l5[. Here the surface atoms are displaced 
in exactly the same way as the elevation maxima of the 
surface profile in the present study. In terms of 2D space 
groups the transition is from p4 to p2mg implying a dou- 
bling of the unit cell. The rectangular p2mg symmetry 
calls for a rectangular metric, that is for different lattice 
parameters along two orthogonal directions but this has 
not observed in our experiment, presumably because of 
the insufficient resolution. 
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FIG. 1. Neutral stability diagram for an excitation fre- 
quency Q = 27r X 9.5 Hz slightly above the bicritical point 
— 2ti X 8.7Hz. Within the left (right) resonance tongue 
plane wave perturbations (wavenumber k) with a subhar- 
monic (harmonic) time dependence become unstable. The 
absolute minimum (ac — 4.25^, ks — 423m~"'^) determines 
the threshold for the onset of Faraday waves. The dashed 
horizontal lines indicate the drive amplitudes at which the 
patterns of Figs. 3 and 4 are taken. 
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FIG. 2. Phase diagram of the observed nonlinear pat- 
terns. Region /: fiat surface- //: subharmonically oscillating 
squares as depicted in Fig. |3|a,b. ///: ^2 x>/2 superlattice 
(Fig. ^,b); IV: VS Xa/3 superlattice (Fig. ^,f); V: harmoni- 
cally oscillating hexagons (Fig. ^,h). The symbols mark the 
experimental observed transition points between the phases. 
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FIG. 3. a,b): Photograph and Fourier spectrum of the 
subharmonic square pattern (region II in Fig. |l|). All video 
images show the central region of the cylindrical container, 
whose diameter is about twice as large, c) video image com- 
puted by the ray tracing technique d) density plot of the as- 
sociated 2D surface profile. 
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FIG. 4. Photographs and related Fourier spectra taken 
at drive amphtudes as indicated in Fig. 1; a,b): y/2 x a/2 
superlattice corresponding to region III of Fig. jij; c,d): 
"quasi- hexagonal" transient at the cross-over between /// 
and IV; e,f): y/S Xa/3 superlattice (region IV); g,h) hexag- 
onal pattern (region V). 



